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Abstract 

A large volume TPC is a leading candidate for the central tracking detector at a 
future high energy linear collider. To improve the resolution a new readout based on 
micro-pattern gas detectors is being developed. Measurements of the spatial resolu- 
tion of cosmic-ray tracks in a GEM TPC are presented. We find that the resolution 
suffers if the readout pads are too wide with respect to the charge distribution at 
the readout plane due to insufficient charge sharing. For narrow pads of 2 x 6 mm^ 
we measure a resolution of 100 /um at short drift distances in the absence of an 
axial magnetic field. The dependence of the spatial resolution as a function of drift 
distance allows the determination of the underlying electron statistics. Our results 
show that the present technique uses about half the statistical power available from 
the number of primary electrons. The track angle effect is observed as expected. 

Key words: Gaseous Detectors, Position-Sensitive Detectors, Micro-Pattern Gas 
Detectors, Gas Electron Multiplier 
PACS: 29.40.Cs, 29.40.Gx 



1 Introduction 



The time projection chamber (TPC) [1,2] has been a mainstay of large particle 
detectors since its initial concept in the 1970's. The traditional TPC has an 
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end cap detector that uses anode wires for amplification of the signal. When 
operated in an axial magnetic field, this leads to the so called E X B effect [3] 
close to the wires, which significantly degrades the resolution of the TPC. 
Proposals to readout TPC signals without the use of anode wires have been 
suggested [4, 5] in the past. The recent development and success of micro 
pattern gas detectors (MPGD) such as the //Megas [6] and the GEM [7,8] has 
renewed interest in this solution. 

The advantages of MPGD detectors are that they require less mass for con- 
struction, should not have any E X B effect, naturally suppress positive ion 
feedback and allow more freedom in the shape and orientation of the readout 
pads. In addition the signals are faster, potentially improving the double track 
resolution in drift time. In the case of MPGDs, the signal on the readout pads 
can be a direct electron collection signal or an induced signal. The advan- 
tage of direct signals is that their amplitude is larger and the signal is more 
confined, thus potentially improving the spatial double track resolution. The 
disadvantage of the confined signal is that the pads need to be much narrower, 
on the order of the width of the ionization charge distribution, which increases 
the number of channels and thus the cost. In the case of GEMs the ionization 
charge can be spread naturally in the gaps between the GEMs and readout 
pads. It is also possible to use the induced signal [9, 10] which has a wider 
spread than the direct signal, but a reduced amplitude. 

GEM amplification with pad type readout planes has been shown to give 
excellent spatial resolution for point sources, such as X-rays converting in 
a gas [11], which is useful for medical applications, where the pad size can 
be arbitrarily small to give the required resolution. In the case of a large 
scale experiment using a TPC, such as the proposed TESLA detector, the 
pad size determines the number of channels and thus the cost; in that case 
it is important to make the pad size as large as possible consistent with the 
resolution required. 

In earlier studies [10], using a double GEM amplification stage, we determined 
the point resolution, s, that can be achieved for X-rays converting in the gas us- 
ing the direct charge distribution near the edge of hexagonal pads (s ~ 70 /xm) 
and the induced charge distribution near the middle of pads (s ~ 80 fim). A 
subsequent study [12] with cosmic rays and a small TPC with an end cap 
detector with 5 staggered rows of 2.5 x 5 mm^ rectangular pads showed that 
these pads produced an adequate track resolution using the direct charge. 

In this paper we examine the resolution that can be achieved using the direct 
signal from a double GEM amplification stage and a rectangular staggered 
pad readout scheme. In particular we examine the effect of the pad width and 
length on the spatial resolution and attempt to gain a better understanding of 
the phenomena that affect the resolution. For this purpose we measured the 
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spatial resolution as a function of several different quantities, including three 
different pad sizes and local position across a pad, two gases, drift distance, 
crossing angle, and signal amplitude. 

The two gases used were PIO (Ar(90):CH4(10)), a fast gas with large diffusion, 
and Ar(90): 002(10), a slow gas with relatively small diffusion, operated at a 
voltage below the peak velocity. The different diffusion properties allowed 
us to study the effect of pad size relative to the width of the direct charge 
distribution arriving at the pads, and to simulate, with the Ar002 mixture, 
reduced diffusion as would be present with a PIO type gas in a magnetic field. 



2 Experimental setup 

The test TPC used for these measurements is housed in a cyhndrical pressure 
vessel filled with PIO or Ar002 gas at atmospheric pressure. The TPO has 
a maximum drift length of 15 cm and an active area of 8 x 8 cm^. The drift 
field of 138 V/cm is established by a series of thin window frame electrodes 
located between the cathode plane at the far end and the readout end plane 
at the other end of the TPO. A charged particle crossing the drift region will 
ionize the gas; the released electrons drift to the end plane where they are 
amplified and detected on a readout POB. While drifting the charge cloud 
gets wider due to transverse diffusion, an effect that would be reduced in an 
axial magnetic field. 

We use a double GEM structure for amplification with a gap of 2.4 mm be- 
tween the first and the second GEM. The voltage difference across this transfer 
gap is 653 V resulting in a field of 2.7 kV/cm. The induction gap between the 
second GEM and the readout board is 5.4 mm wide with a voltage difference 
of 1783 V and a field of 3.3 kV/cm. The transfer and the induction gaps were 
purposely large to diffuse the electron cloud and thus spread the signal over 
more than one readout pad. 

The OEMs were made from 50 //m thick kapton foil coated with copper on 
both sides. The holes with a diameter of ~90 //m at the surface are arranged 
in a hexagonal pattern with a spacing of dhex = 140 /im. The voltage across 
the OEMs is 357 V each. Gain measurements were made for 5.9 keV ^^Fe x-ray 
conversion electrons in the gas and used standard pulser calibration technique. 
The effective gains for Ar002 and PIO were about 6700 and 4600 respectively. 

The readout-pad layout shown in Figure 1 consists of 192 pads which are 
reduced to 64 readout channels via multiplexing. The pads in the outer rows 
(1,2,4,7,9,10) are 2.54 mm x 5.08 mm large; in the test row 5 the pads are 
2.032 mm x 6.096 mm large and in row 6 they are 3.048 mm x 5.080 mm 
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large. Rows 3 and 8 consist of wide pads used for filtering. The outer pads in 
rows 1-4 and 7-10 arc multiplexed to give one veto channel on the left and 
right side, respectively. 

We use a right-handed coordinate system with the x-coordinate horizontal 
and the y-coordinate pointing upwards; the z-coordinate corresponds to the 
drift distance with 2; = at the first GEM. The azimuthal angle (f) and the 
polar angle 9 are measured with respect to the y-axis. 

The signals are read out via ALEPH TPC wire preamplifiers [13] and 64 
channels of 200 MHz, 8 bit EADCs custom made at the University of Montreal. 
For data acquisition we use the MIDAS [14] system. 

A three layer scintillation counter telescope is used to trigger on cosmic-ray 
tracks. One scintillator counter is placed above the TPC and two below, sep- 
arated by a 10 cm thick layer of lead. The ~19 cm width of the counters and 
the distance of ~40 cm between the two outer counters defines the acceptance 
coverage in z. 



Run 1010, Event 22 




I 

-30 -20 -10 10 20 30 4 

X [mm] 



Fig. 1. Event display with a reconstructed track in ArC02 gas. The shade of the 

pad corresponds the the reconstructed signal ampHtude. The lower threshold is 2 
and one hit per row of more than 8.5 is required for the track fit. In general there 
is a 3-fold multiplexing for the outer and test rows. The rows are numbered from 
bottom to top. 
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3 Theory 



The observed width ^ of the track, (Jtrack, which is the standard deviation of 
the charge cloud perpendicular to the track, is determined by two components, 
the spread associated with the readout system ctq and the transverse diffusion 
depending on the drift distance z: 



C^track = O-0 + C'd ^ (1) 
2 _ 2 I 2 I 2 

•^0 ~ "^hex "T '^intern "r '-'other ) 

where Cd = y^lDxjv is the 1-dimensional diffusion coefficient given by the 
transverse diffusion constant and the drift velocity v. In a magnetic field 
Dt{B) = Dt{0)/{1 + cj^ r^), thus resulting in a reduced transverse diffusion. 
The contribution (Jq is composed of several parts. The first term originates from 
the hexagonal pattern structure of the GEM depending on the hole distance 
dhex- For our geometry (Jhex is estimated to be ~ 50 /im. The second term, 
Cintern, results from diffusion between the GEMs and the readout pads. For 
the present setup o-mtcrn — 318 fim for ArC02 and ~ 460 //m for PIO. Other 
factors denoted by (Tother also contribute. 

The standard deviation cr^ of the charge cloud distribution across a row of 
pads also includes the crossing angle effect o",^ depending on the track angle 
and is given by: 

^l = ^l + (^l + ^l (2) 
(^D — C'd cos 
a^ — L/Vvi tan0 . 

The factor 1/ cos0 in the transverse diffusion term ctd reflects the projection 
of the charge distribution onto the x-axis. The crossing angle effect comes from 
the spread of xiy) for a track with an angle (p. Projected onto the x axis this 
leads to a rectangular function of total width L tan 0, where L is the length of 
the pad. The standard deviation for such a rectangular uniform distribution 
is given by a^. 

The observed x-resolution Sx reflects the precision with which the mean of the 
charge distribution can be determined and hence has additional factors from 
statistics: 



^ Throughout this paper the width of a distribution refers to its standard deviation. 
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2 _ 2 I 2 I 2 

*0 ~ '^hex + '^intern ^other 



Most contributions depend on the number of electrons Ut produced by the ion- 
izing particle. Some of these electrons stem from secondary ionization. They 
are therefore correlated forming rid clusters. The number of electrons and clus- 
ters created across a row of pads is Nt — rif Lj cos0 and N^x = rid ■ L/ cos^, 
respectively. For example for Argon rid — 24.3/cm and = 94/cm [15]. 

All components of sq depend on the signal amplitude. The contribution from 
the GEM structure Shex is minor: the track width due to diffusion is large 
enough to cover a sufficient number of GEM holes. The effect from internal 
diffusion Sintem is strongly reduced due to the high gain. Contributions from 
electronic noise, calibration errors and limitations due to insufficient charge 
sharing between the pads are included in Sother- 

The effect from transverse diffusion depends on the effective number of elec- 
trons N^^ — R- N^^, which includes a reduction factor R. The crossing angle 
effect depends on the effective number of clusters N^^ = {NdY. Following the 
notation of [16] the number of clusters is reduced by the exponent e. 
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Fig. 2. Determination of time and amplitude of a pulse, see text. The dot indicates 
the reconstructed TO and amplitude. 
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4 Reconstruction 



The analysis package is based on Fortran £95 code [17]. In a calibration run 
pedestals and pulse fall times ifaii as well as the relative gain are determined 
for each readout channel. The time TO and the amplitude of the signals are 
determined from the pedestal subtracted ADC pulse as shown in Figure 2. The 
time Tpeak is determined as the time bin with minimum ADC count. In the re- 
gion [Tpeak + 50; Tpeak + 350] an exponential A{t) = Apeak *exp -{t - Tpeak) Afaii 
is fit to the ADC spectrum to determine the amplitude Apeak at Tpeak- The 
time TO is determined via a linear fit to 25 time bins at the rising edge as 
ADC(TO) = Apcak/2 and the signal amplitude is —A{TO). The amplitudes are 
corrected for the relative gain of each channel. The RMS of the correction 
coefficients is 5%. Only signals with an amplitude of more than 2 ADC counts 
are recognized as pad hits and taken into account. Events are rejected if a veto 
channel has an amplitude of more than 8.5 ADC counts. The TO of a row is 
determined as the amplitude weighted mean of the times of the hits in this 
row. The sum of reconstructed amplitudes in row 5 is shown in Figure 3. 

The track fit is performed similar to [18]. In the upper and lower two rows (1,2 
and 9,10) start points arc determined from a centroid calculation of the largest 
amplitude channel and its neighbor pads. These two points arc connected by a 
line to form the seed track. Because of the multiplexing several seed tracks are 




amplitude 



Fig. 3. Sum of pad signal amplitudes in row 5 for both gases. The distribution 

is proportional to the energy loss, the difference between ArC02 and PIO is due 
to different gain. The RMS of the noise on a single pad is ~0.5 ADC counts with 
typically 2 to 3 pads contributing to the amplitude of the row. The cutoff for a pad 
signal was 2 ADC counts. 
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Table 1 

MAGBOLTZ (version 5.1) predictions for a drift field of 140 V/cm. 





drift 


diffusion 


gas 


velocity 


transverse 


longitudinal 




(//m/ii.s) 


(//m / v/cm) 


(//m / ^/cm) 


PIO 


55 


564 


374 


ArC02 


8.9 


229 


241 



found and the track with the most rows having a related hit with an amplitude 
of more than 8.5 is chosen. In general this choice is unique. There must be at 
least six rows with hits out of the eight outer and test rows. Other events are 
rejected from the analysis. 

The track projection in the x-y plane can be described with three parameters: 
the x-position at y=0, xo, the track angle, and the width of the charge 
cloud, (Ttrack- Thc track parameters are determined from a maximum likeli- 
hood fit where a uniform line of charge with a Gaussian profile is assumed. 
This idealized distribution is integrated over the pad areas and normalized 
across a row to obtain the expected charge probabihties. Prom these and the 
observed signal amplitudes a likelihood function is determined, which includes 
a uniform noise probability of 0.2%. The noise level is determined from the 
data; a variation between 0.1 and 0.5% has only a small effect on the fitted 
track parameters. 

The drift distance at y=0, zq, and the angle 6 are determined from a straight 
line fit to the TO of each row as a function of y. All eight rows are used to 
determine global distributions of track angles (p and 9 as well as Xq and Zq. 
The drift velocities as determined from the data are 55±4 //m/ns for PIO and 
8.3±0.3 //m/ns for ArC02. The result for PIO is in good agreement with the 
prediction from MAGBOLTZ given in Table 1, while the measured velocity for 
ArC02 is smaller than expected. This might be due to a limited time window 
recorded or because of gas impurities. In ArG02 we lose about 10% of the 
electrons over the full drift distance of 15 cm due to attachment. This is an 
indication of impurities in the gas which may affect gas properties. No such 
effect is observed with PIO. 



5 Analysis Results 

The resolution for the drift distance z is shown in Figure 4 as a function of 
drift distance. The intrinsic time resolution is about 13 ns for PIO and 9 ns 
for ArC02. It is worse for PIO since the average signal amplitude is smaller. 
While the z resolution for PIO is completely dominated by the intrinsic time 
resolution, the effect of longitudinal diffusion is visible for ArC02 because of 
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Fig. 4. Resolution of drift distance for both gases as a function of drift distance for 
small track angle \(/)\ < 5°. 

the much smaller drift velocity. The observed dependence is linear and not 
quadratic as expected: s^/iiTa = 80 + 14 * z/cm. It does not depend on the 
readout pad size. 

The width of the charge cloud Utrack is shown in Figure 5 as a function of the 
drift distance. The mean transverse diffusion coefficient can be determined 
from a hnear fit to (jfrack(-^) (Equation 1). We obtain 

Cd = 0.429 mm/ Vcm for PIO and (4) 
Cd = 0.209 mm/^/aR for ArCOa , (5) 

with negligible statistical errors. The result for PIO is smaller than the expec- 
tation from MAGBOLTZ given in Table 1 while the result for ArC02 is in 
reasonable agreement. The width of the charge cloud at z = is determined 
to be 0.563 ± 0.006 mm for PIO and 0.544 ± 0.006 mm for ArCOs. Only a 
part of this measured amount is expected from o"intcrn- Since o"hcx is negligible 
a contribution of several 100 /^m remains unexplained for both gases and must 
be assigned to other factors in Uother- 

The remainder of this paper concentrates on the study of the x resolution 
dependence on track angle, transverse diffusion and amplitude. In this analysis 
the track parameters are not determined from reconstructed hit-positions in 
each row but from a fit to the charge distribution of the full track. Therefore 
the concept of the point position in a row is not a priori given. The x-position 
in a row, Xj-ow is determined from a track fit to the charge distribution in this 
row only, keeping all track parameters fixed apart from x. The x-resolution 
Sx is obtained from the width cr of a Gaussian fit to the distribution of the 
residuals 5 — Xrow — s^track; Xtvack — xq + tan * 7/row K these residuals are 
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Fig. 5. Fitted track width Oxxz-cV as a function of drift distance z. The upper plots 
show CTtrack versus z for all events, the lower plots show the average cr^^ack versus z. 

derived from a track fit including the test row the obtained spread o"" will be 
systematically too small. On the other hand the spread cx'^^ obtained from a 
track fit excluding the test row will be too large. As shown in appendix A the 
geometric mean of these two results — -\/a^°*~a^ gives the correct estimate 
for the point resolution. For these studies the charge width Utrack is fixed to 
the mean observed track width as a function of the drift distance. 



5.1 X resolution depending on pad width 



First we investigate the dependence of the resolution on the width of the 
pad. To eliminate other effects only tracks with |0| < 5° are used for this 
study. Figure 6 shows the x resolution in 3 mm wide pads as a function of 
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Table 2 

Mean track width and RMS for three regions of drift distance and two gases. 





track width (mm) 


drift distance 


ArC02 PIO 


- 3 cm 


0.53±0.14 0.72±0.24 


3-8 cm 


0.67±0.13 1.12±0.21 


8-15 cm 


0.81±0.11 1.52±0.21 



the distance between the reconstructed position and the centre of the pad 
for ArC02. To obtain samples with different diffusion, i.e. size of the charge 
cloud, three regions of drift distance are considered as given in Table 2. For 
short drift distances, hence small charge-cloud size, the resolution gets signif- 
icantly worse in the centre of the pad. This is because an increased fraction 
of signals is collected only on one pad and charge sharing is less effective for 
the determination of the position of the track in this row. This effect leads 
also to a non-uniform distribution of the measured Xy-ow positions, where more 
hits are reconstructed in the center of a pad if the pad is too wide. This study 
is repeated for the 2 mm and 2.5 mm (row 4) wide pads, which shows that 
the effect sets in if the pad is wider than about three times the width of the 
charge cloud and becomes prominent for a pad width larger than four times 
the charge width. However, this effect depends also on the amplitude; signals 
with large amplitude tend to have more charge sharing. All measurements are 
made without magnetic field. However, the results can be reinterpreted for a 
given width of the charge cloud. This study indicates that 3 mm wide pads 
are too wide for a charge cloud with a width of less than about 1 mm, which 
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Fig. 6. X resolution as a function of distance to the centre of the pad for ArC02 
and three regions of drift distance. The points are results from 3 mm wide pads 
with |(^| < 5°. 
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is the case for ArC02 and PIO at small drift distances. To avoid this effect, 
the following analyses are restricted to the 2 mm wide pads. 



5.2 X resolution depending on drift distance z 



The resolution deteriorates with increasing drift distance because of transverse 
diffusion. This is studied using tracks with small angle, \4>\ < 5°, to suppress 
the track angle effect. As can be seen in Figure 7 the effect is less pronounced 
in ArC02 because of the very small diffusion. The function 



2 I 



(6) 



can be fit to the data. Using Cd as given in Equations 4 and 5 we obtain Nf = 
19±7 for ArC02 and 20.6±0.7 for PIO. These numbers are smaller than the 

total expected number of electrons Nt. However, the wide range of amplitudes 
makes the interpretation difficult. The expected A^t as given in [15] relates to 
the mean amplitude observed in a row, while the resolution is proportional to 
1 / \/At , which results in a bias. For the full range the mean amplitude (A) is 
larger than l/( (I/VA) f by 30% for both gases. 
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Fig. 7. Resolution function of drift distance for 2 mm wide pads, 

tracks with < 5°, both gases. 
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If the sample is split up in three regions of signal amplitude this bias is reduced 
to about 10%. The amplitude ranges are chosen such that the number of events 
in each range is approximately the same. The fit results are summarized in 
Table 3 and Figure 8. The ratio R is consistent with 0.5 independent of the 
amplitude. This interpretation does not account for the ~10% loss of primary 
electrons over the full drift-distance observed in ArC02. There is no electron 

Table 3 

Fit result for resolution as a function of drift distance. The number of electrons A^t 
for tracks with small angle \(p\ is scaled with mean amplitude; for the full range rit 
is taken from [15]. The errors are statistical only. 



Amplitude 








derived 


range 


mean 




{fim 1 ^cm) 




^rcff 

t (%) 


PIO 












all 


93 


55 


94.6±1.6 


20.6±0.7 




< ^ < 60 


42 


25 


122.3±2.6 


12.3±0.5 


50.0±2.1 


60 < A < 100 


78 


46 


93.9±2.0 


20.9±0.9 


45.3±1.9 


100 < A < 1000 


175 


103 


70.0±3.1 


37.6±3.3 


36.3±3.2 


ArCOa 












all 


150 


57 


47.8±9.3 


19±7 




< ^ < 100 


69 


26 


55.6±3.0 


14ib2 


54±6 


100 < ^ < 170 


131 


50 


40.2±2.8 


27±4 


55±8 


170 < A < 1000 


279 


106 


23.2±3.9 


82±27 


77±26 



0.8 



0.6 



0.4 



0.2 




0.8 



0.6 



0.4 



0.2 




100 200 300 50 100 150 200 

amplitude amplitude 



Fig. 8. Ratio R = /N^ as determined from the resolution as a function of the 
drift distance for both gases and 3 regions of amplitude. 
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loss in the drift-region for PIO. No loss of GEM transparency is suggested at 
our operating voltages from published analyses [19]. However, it is clear that 
the present resolution measurement technique effectively uses only about half 
the statistical power available from the number of primary electrons. 



5.3 X resolution depending on track angle (f) 



The track angle effect on the resolution is expected to be oc tan and depends 
on the length of the pad as well as the effective number of clusters N^i. Figure 9 
shows the x-resolution as a function of track angle for drift distances of less 
than 3 cm for three regions of amplitude. The following function is fit to the 
data: 



«x = \^' + T7Ttan2(0-(^)/7V^^i, (7) 



where (p is an additional free parameter allowing for a bias in the track angle 
and s includes contributions from diffusion sd- For the fit the number of 
clusters Nd is taken to be independent of the amplitude, all variations being 
assigned to the exponent e. The fit results for s, (p and e are given in Table 4. 

The offset cp is consistent with indicating that there is no systematic shift. 
The number of primary clusters A^ci is reduced to the number of effective 
clusters by the exponent e. We see no dependence of e on the amplitude. As 
expected the basic resolution s improves with amplitude. Due to the high gain 
s is mainly determined by the diffusion Sd and not so much by the internal 
resolution sq- 



Table 4 

Fit result for track angle effect. 



Amplitude 




Fit result 




range 


s (mm) 


f (deg) 


e 


PIO 








< ^ < 60 


0.170±0.005 


-0.3±0.5 


0.50±0.03 


60<A< 100 


0.122±0.004 


-0.5±0.3 


0.54±0.02 


100 < A < 1000 


0.111±0.004 


0.3±0.3 


0.49±0.02 


ArC02 








< ^ < 100 


0.130±0.005 


0.1±0.5 


0.64±0.03 


100 < ^ < 170 


0.103±0.003 


0.2±0.3 


0.56±0.02 


170 < ^ < 1000 


0.105±0.003 


0.1±0.3 


0.55±0.02 
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Fig. 9. X resolution as a function of track angle for 2 mm wide pads and drift 
distance less than 3 cm. For both gases the sample is split up in three regions of 
amplitude A. The expected angular dependence is fit to the distributions. The fit 
results are given in Table 4. 
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6 Conclusion 

We have investigated the dependence of the spatial resolution of cosmic-ray 
tracks in a GEM TPC on various parameters. We found that the resolution 
degrades if the pads are too wide with respect to the track-charge width ar- 
riving at the readout plane due to insufficient charge sharing between readout 
pads. The observed charge width is larger than that expected from transverse 
diffusion between the GEMs and the readout plane. Therefore 2 mm wide 
pads are large enough to achieve a resolution of 100 /xm for drift distances of 
less than 3 cm in the absence of an axial magnetic field. The dependence of 
the spatial resolution as a function of drift distance allows the determination 
of the underlying electron statistics. Our results show that with the present 
technique effectively only about half of the primary electron statistics is used 
for the position determination. The track angle effect is observed as expected. 
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A Corrections to the x resolution 



When determining the point resolution the typical method involves fitting 
a straight line to all points and determining the standard deviation of the 
residuals; this method gives a resolution that is too good, since the point for 
which the resolution is to be determined was included in the line fit. The 
alternate method is to fit a straight fine without the point for which the 
resolution is to be determined; this gives a resolution which is worse than 
the actual resolution since the line is determined from the other points which 
themselves have an uncertainty. 

A detailed analysis reveals that a better estimate of the true resolution is given 
by the geometric mean of the two methods, that is: af — cr^m • cr^ex, where CTj is 
the better estimate of the resolution for point i, 5™ and Sf^ are the measured 
residuals when the point is included and excluded respectively, and cr^in and 
(T^cx are the standard deviations of the residual distribution when the point is 
included and excluded from the fit. 

Let us assume a track consisting of N measurements with known values yj, 
1 < J < The corresponding measured values Xj are distributed around the 
expected mean (xj) = a + b yj with the standard deviations sj, where a and 
b are the track parameters. To determine the resolution of one measurement 
i it is convenient to choose the coordinate system so that yi = 0. In this 
case, the residual is given by 5i — a — Xi, where a can be determined from 
a least square fit to the track by either including (a™) or excluding (af^) the 
measurement i. The residual 6i will be distributed with a standard deviation 
(T5. which is related to Sj, but depends on the coordinates {xj,yj) and weights 
Wj = of all measurements. 

Minimising the gives an estimate for a: 



a- 



m 



Ej WjXj ■ Efc Wkyl - Efc Wfci/fc ■ Ej WjXjyj 



, where 



(A.1) 



D 



iin 



and 



ex 



Ej^i wjxj ■ Ek^i Wkyl - Ek^i WkVk ■ Ej^i WjXjyj 



, where 




lex 
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And since y{j=i) — 

Dr = D^''-w,J2^,y]. (A.3) 



The residual 5™ of point i is: 
5]-^ = a'" - Xi 



Ej^i WjXj (Efc Wkvl - Vj Efc WkVk) - DT^i 



Assuming that the measurements are independent, the variance of the resid- 
ual distribution is approximately: 



'dSl 



4, = E(^^.l ■ (A.5) 



The partial differentiation picks out the Xj terms yielding: 



Ej^i w s (Eit Wkvl - Vj Efc WkVk) - {DfSi) 



Expanding, rearranging and collecting terms yields: 

(T|n = s- and similarly: ct|x = s,^— . (A.7) 

The quantities and D^^ are fixed for a given layout and can be calculated 
to correct the resolution measured, however it is simpler to combine the last 
two expressions: 

— agin • cr^ex . (A. 8) 

With this one can thus get a better estimate of the resolution by taking the 
geometric mean of the resolution as determined with the point included in and 
with the point excluded from the fit without having to calculate a correction 
factor. As expected, our tests show that for a large sample the results are 
identical for the two methods. 
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